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Abstract 



We study the asymptotic behavior of the ghost propagator in the quenched 
Q-i, SU(3) lattice gauge theory with Wilson action. The study is performed on 

lattices with a physical volume fixed around 1.6 fm and different lattice 
spacings: 0.100 fm, 0.070 fm and 0.055 fm. We implement an efficient 
algorithm for computing the Faddeev-Popov operator on the lattice. We 
are able to extrapolate the lattice data for the ghost propagator towards 
the continuum and to show that the extrapolated data on each lattice 
can be described up to four-loop perturbation theory from 2.0 GeV to 6.0 
GeV. The three-loop values are consistent with those extracted from pre- 
vious perturbative studies of the gluon propagator. However the effective 
A.Ms" scale which reproduces the data does depend strongly upon the or- 
der of perturbation theory and on the renormalization scheme used in the 
parametrization. We show how the truncation of the perturbative series 
can account for the magnitude of the dependency in this energy range. The 
contribution of non-perturbative corrections will be discussed elsewhere. 
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1 Introduction 



Whereas lattice gauge theory (LGT) has been initially formulated to study gauge- 
invariant quantities in the non-perturbative regime, it has long been recognized 
that LGT could be a useful tool for studying gauge-variant quantities such as 
Green functions, both in the non-perturbative and in the perturbative regimes. 
The SU(3) gluon propagator in momentum space was first considered P to gain 
some insight into the physics of confinement. Much work was then devoted to 
the study of its infrared behavior (for a review see PJ. Subsequent studies 
[3J m were focused on the ultraviolet behavior and have been able to compare 
quantitatively the large momentum dependence of the lattice gluon propagator 
with perturbative predictions beyond one-loop order. The result for A^jg was 
found to depend strongly upon the order of the perturbation theory and upon 
the renormalisation scheme used in the parametrization. This strong dependence 
raised the question whether the energy windows in these calculations were large 
enough for perturbative QCD to be a valid approximation. 

On the other hand, as shown by Gribov [3] , the infrared behavior of the gluon 
propagator is closely related to the singularity structure of the ghost propagator 
inferred from the gauge-fixing ambiguities. As is well-known, the Landau gauge, 
which is presently the only covariant gauge for which there exists effective lo- 
cal gauge-fixing algorithms on the lattice, suffers from these ambiguities. The 
comprehensive theoretical study by Zwanziger |]] of the Faddeev-Popov operator 
on the lattice in Landau gauge spurred the first numerical study of the ghost 
propagator 7j in SU(2) and SU(3) gauge theories. Most subsequent activity has 
been dedicated to the SU(2) lattice gauge theory in the infrared region, mainly 
for technical reasons as we shall explain below. There are relatively few numer- 
ical studies of the SU(3) ghost propagator which are either more focused on the 
infrared region and the Gribov copy problem [HJ ED] or have only performed a 
qualitative perturbative description in the quenched approximation fTTJ H2] and, 
quite recently, in the unquenched case also [T3] . 

It is important to make the study of the SU(3) ghost propagator in the ultravi- 
olet region more quantitative for comparison purposes with the gluon propagator. 
Lattice results at small distances may be described by perturbation theory and 
the independent extraction of the Aq CD scale from the two propagators would 
provide a self-consistency test of the analysis and of the lattice approach. It 
would be particularly significant to confirm or not, from the study of the lattice 
propagators alone, the need for the non-perturbative power corrections found in 
the study of the three-gluon coupling on the lattice [Ti] . 

The paper is organized as follows. We will begin by recalling in section El 
the method used to relate lattice data for the ghost propagator to its pertur- 
bative renormalization description. Then we proceed by exhibiting in section |3] 
the salient features of our lattice calculation, particularly of our implementation 
of the Faddeev-Popov operator on the lattice. The following section outlines the 
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general method that we devised previously Hj to eliminate hypercubic artifacts 
from two-point functions and extrapolate the lattice data towards the continuum. 
This extrapolation is crucial to succeed in a quantitative description. The results 
are discussed in section which contains several subsections where the analysis 
is performed in different renormalization schemes up to four-loop order. In par- 
ticular the scheme dependence is thoroughly investigated and used to probe the 
effects of the truncation of the perturbative series. We conclude in section El with 
a comparison of the different methods to compute the Aqcd scale on the lattice. 



2 Renormalization description of the ghost prop- 
agator 

Let be some gauge-fixed multiplicatively renormalizable one-particle irre- 
ducible n-point bare Green function defined in euclidean momentum space and 
in some regularization scheme with cut-off A. Let s denotes some polarization 
state and kinematical configuration of the external particles contributing to Tg. 
Let p denote a scale transformation on s and gs denote the bare coupling. It is 
well known that, in any renormalization scheme R defined by some renormaliza- 
tion conditions on state s at the renormalization point p = p, we have 

ri n) (p, s, g B , A) = Z TtR (ji, s, g R , A)r£>(p, s, g R , p) + ^A" 1 ) (1) 

where Z^ >R is the renormalization constant in scheme R } is the renormalized 
Green function and g R {p) is the renormalized coupling. We omit the dependence 
on the gauge parameter for simplicity of notation since we will specialize to 
Landau gauge. 

The explicit dependence on p drops out of the renormalized Green function 
at the renormalization point p = p. It follows that 

^ rfln (rjW^ = ^ dln(Z T , R (p,s,g R ,A)) + dln(T%(s, g R )) 



A^oo d In p 2 A-^oo d In p 2 d\np 

= Jr,n(9R) + 



dg R d\nT n R 



d\np 2 dg R 

(2) 

The arbitrariness in the choice of the renormalization scheme R has prompted 
attempts at determining the "best" schemes for describing the g 2 -evolution of 
bare Green functions on the lattice. Clearly it is always possible to find a change 
of coupling which will be a best approximation of a set of data at a given order of 
perturbation theory, within some prescribed criteria. Rather than pursuing this 
route, we will follow the standard wisdom which consists in choosing renormal- 
ization conditions appropriate to the continuum quantity under scrutiny. 
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Momentum substraction schemes have long been used to define renormaliza- 
tion conditions befitted to the description of the renormalization dependence of 
"physical" quantities. They are defined by setting some of the 2 and 3-point 
functions to their tree values. In the MOM schemes, for these Green functions, 
Eq. (J2J) simplifies to 

r dhi (ltV, g ,^,A)) dHz m) 

Jim -r ^ = -r 2 = 7r,MOM{9MOM) (3) 

a^oo d In \i l a In \i A 

Infinitely many MOM schemes can be defined which differ by the substraction 

point of the vertices. We have shown in [Hj that the MOMg scheme defined by 
substracting the transversal part of the three-gluon vertex at the asymmetric 
point where one external momentum vanishes, appears to provide a much better 
estimate of the asymptotic behavior of the gluon propagator in Landau gauge 
than the MS scheme. For the study of the asymptotic behavior of the ghost 
propagator in Landau gauge, it seems therefore natural to use a MOMc scheme 
defined by substracting the ghost-gluon vertex at the asymmetric point where 
the momentum of the external gluon vanishes. Comparison of the two MOM 
schemes should provide us with an estimate of the systematic error entailed in 
the truncation of the perturbation theory. 




Figure 1: MOMg scheme (left) and MOMc scheme (right). 



The perturbative calculation of the gluon, ghost and quark self-energies and 
all 3-vertices appearing in the QCD Lagrangian have been done at three-loop 
order in the MS scheme and in a general covariant gauge at the asymmetric point 
with one vanishing momentum [T^]. These three- loop results allow one to relate 

the coupling constants of any MOM-like scheme to the MS scheme at three-loop 
order. For the MOMg and MOMc schemes defined above these relations read 
respectively in Landau gauge and in the quenched approximation (n/ = 0), with 
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16vr 2 



i,„ 70 h * , ^16217 _ 153 \ 3 

ft MOM g - "ms + 3 Ms + I 576 4 ^ i M 

/304676635 299961 81825 \ 
+ I - 6912 64~ Cs " ~6T C6 J * 



(4) 



223 2 / 918819 _ 351 \ 3 
Vomc ~~ MS + + ( 12 96 8~ J m 

/ 29551181 _ 137199 _ 74295 \ 4 
+ {~ 864 32" C3-_ 64 _Cs) ^ 



(5) 



MS 



The very large coefficients of these perturbative expansions explain the difficulties 
met by the MS scheme to approach asymptotic scaling below 10 GeV. 

The recent calculation [T£| of the anomalous dimensions in the MS scheme of 
the gluon and ghost fields at four-loop order, together with the knowledge of the 
/3-function |17j . makes it possible to perform the analysis of the lattice data for 

the gluon and ghost propagators up to four-loop order also in the MOMg and 
MOMc schemes. The numerical coefficients of the /5-function defined as 

P(h) = -^- 2 =-J2 &h i+2 + C(/i n+3 ) (6) 



dlnu 2 

r i=0 



arc: 



^MOMg = 24m6) ^momc = 2952.73, /3 3 MOMg = 84353.8, #f OMc = 101484. (7) 

For completeness we also give the expansion coefficients of the renormalisation 
constants of the gluon and ghost fields in the MOM schemes with respect to the 
renormalized coupling of the MS scheme up to four-loop order: 



rfln/x 2 ~ 2 Vs+ 24 ^ + V 



,3 



d\n{Z^MO M ) _ 13 , 3727, 2 /2127823 9747 ..., 
h m + ~TTT~ % + I — ggg ^ps j n- MS 

( 3011547563 18987543 1431945 . j i 
+ V 6912 206 - Cs 6T~ C 7 Ms 

d\n(Z 3MOM ) _ 9 , , 813, a /157303 5697. \ 3 
ft MS + -77T % + — 7T A ^-<-,3 ri 



(8) 



dln/x 2 4 -ms^ 16 "m~r y 64 32 s:> / - ^ 

/219384137 9207729 221535 , , , 
+ V 1536 512~ C3 32" Cs J MS 



(9) 



,4 
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3 Lattice calculation 



3.1 Faddeev-Popov operator on the lattice 

The ghost propagator is defined on the lattice as 

G{x - y)S* = ((M- 1 )*) (10) 

where the action of the Faddeev-Popov operator M on an arbitrary element rj 
of the Lie algebra SU(N) of the gauge group SU(N), in a Landau gauge fixed 
configuration, is given by [Hj: 

(M V ) a (x) = 1 J2{ G t& Oft* + P)- v\x)) -(*<-> s - A) 

+ \f abc {r?{x + fiA%{x) - r,\x - fa)A%x - £)) } (11) 

and where, with antihermitian generators T a , 

GfW = \jr{{T\T b } (U,(x) + Ul(x))) (12) 
A^x) = -Jr (T c (U,(x) - Ul(x))) (13) 

Most lattice implementations of the Faddeev-Popov operator have followed closely 
the component- wise Eqs. (jllHlHj) . But the derivation in 6\ shows that the Faddev- 
Popov operator can also be written as a lattice divergence: 

M(U) = -iv • D(U) (14) 

where the operator D reads 

D^{U)r]{x) = i {U^(x)t](x + fa) - i](x)U^(x) + rj{x + /i)C^(x) - Ul(x)r)(x)) 

(15) 

Using conversion routines between the Lie algebra and the Lie group, eqs. (fTH 
IT5j) allow for a very efficient lattice implementation, sketched in Table Q which 
is based on the fast routines coding the group multiplication law. 

3.2 Inversion of the Faddeev-Popov operator 

Constant fields are zero modes of the Faddeev-Popov operator. This operator 
can be inverted only in the vector subspace K 1 - orthogonal to its kernel. If the 
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! Vin, Vout are the ghost fields. 
! U is the gauge configuration, 
type (SUN) [/(*), dU, W, W+, W_ 
type (SUN) r] in (*),r] ou t(*) 
for all x : 

dU = 0. 

W = i] in (x) 

do /i = 1, 4 

W+ = 7] in (x + ft) 

W_ = rj in (x - ft) 

dU = dU + U^x - p) x W + W x U^x) 

- U^x) x W + - W- x U^x - p) 

enddo 

Vout (x) = dU- dU^ - ^MdU - dU^) 



Table 1: Pseudo code of our implementation of the Faddeev-Popov operator. 

Faddeev-Popov operator has no other zero modes than constant fields, then the 
non-zero Fourier modes form a basis of K^: 

r] (x) = Y^ c P eiP ' X > V?7 G K L (16) 

The standard procedure has been to invert the Faddev-Popov operator with one 
non-zero Fourier mode as a source 

S a v [x) = 5 ab e ip - x (17) 

and to take the scalar product of M _1 Sp with the source: 

[S; | M _1 Sp) = (M-X y e- ip - {x ~ v) (18) 

x,y 

= VG(p) (19) 

after averaging over the gauge field configurations. This method requires one 
matrix inversion for each value of the ghost propagator in momentum space. It 
is suitable only when one is interested in a few values of the ghost propagator. 

However, the study of the ultraviolet behavior of the ghost propagator in the 
continuum requires its calculation at many lattice momenta to control the spacing 
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artifacts, as we shall see in the next section. This can be done very economically 
by noting that 

5( X ,y) = l- + ±-Y, e ~ iP<X ~ y) (2°) 



and choosing as source: 



V V 



SZ(x)=5 ab (5(x,0)-±-) (21) 



V 

The Fourier transform of M~ 1 S$, averaged over the gauge configurations, yields 

5>-^(m-^) = Y, e ' ip ' x (( M ' 1 )Z) - ^E e "^(( M_1 ) 



V ^— ' \ v ' x v 



= G(p) - 5(p)J2G(x) (22) 

X 

as a consequence of the translation invariance of the ghost propagator. Therefore, 
with this choice of source, only one matrix inversion followed by one Fourier 
transformation of the solution is required to get the full ghost propagator on the 
lattice. 

There is of course a price to pay, as can be read off Eq. (|22|) which lacks the 
factor V present in Eq. (JTHJ). The statistical accuracy with the source S£ is better, 
especially at high momentum p. However the statistical accuracy with the source 
Sq turns out to be sufficient for our purpose. 

There is one final point we want to make and which has never beeen raised 
to the best of our knowledge. It is mandatory to check, whatever the choice of 
sources, that rounding errors during the inversion do not destroy the condition 
that the solution belongs to K x : 

(M-'S) (x) = (23) 

X 

Indeed, if the zero-mode component of the solution grows beyond some threshold 
during the inversion of the Faddeev-Popov operator on some gauge configura- 
tion, then that component starts to increase exponentially and a sizeable bias 
is produced in other components as well. We have observed this phenomenon 
occasionally, about one gauge configuration every few hundreds, when using the 
implementation of the lattice Faddeev-Popov operator based on Eqs. (|11H13|) . But 
the systematic bias which is induced on the averages over gauge field configura- 
tions can be uncomfortably close to those ascribed to Gribov copies. 

Another virtue of the algorithm described in Table is its numerical stabil- 
ity which is improved by several orders of magnitude. We have never observed 
sizeable deviations from Eq. ()23|) with this algorithm. 



S 



p 


V 


a' 1 (GeV) 


# Configurations 


6.0 


16 4 


1.96 


1000 


6.0 


24 4 


1.96 


500 


6.2 


24 4 


2.75 


500 


6.4 


32 4 


3.66 


250 



Table 2: Run parameters. The lattice spacings are taken from Table 3 in JH] 
with a physical unit normalized by y/a = 445 MeV. 

3.3 The simulation 

We ran simulations of the SU(3) lattice gauge theory with the Wilson action 
in the quenched approximation on several hypercubic lattices, whose parameters 
are summarized in Table 121 All lattices have roughly the same physical volume 
except the 24 4 lattice at f3 — 6.0 which has been included to check out finite- 
volume effects. The SU(3) gauge configurations were generated using a hybrid 
algorithm of Cabibbo-Marinari heatbath and Creutz overrelaxation steps. 10000 
lattice updates were discarded for thermalization and the configurations were 
analyzed every 100/200/500 sweeps on the 16 4 /24 4 /32 4 lattices. 

Landau gauge fixing was carried out by minimizing the functional 

Fu[g] = Re £ £ (l - ±g(x)U»(x)g\x + (24) 

by use of a standard overrelaxation algorithm driving the gauge configuration 
to a local minimum of Fu[g\. We did not try to reach the fundamental modular 
region A, defined as the set of absolute minima of Fu[g] on all gauge orbits. Indeed 
there have been numerous studies, in SU(2) (131201 and in SU(3) |E1IH], of the 
effect of Gribov copies on the ghost propagator. The consensus is that noticeable 
systematic errors, beyond statistical errors, are only found for the smallest p 2 , 
much smaller than the squared momenta that we used to study the asymptotic 
behavior of the ghost propagator. 

Then the ghost propagator G(p) is extracted from Eq. (}2~2~|) for all p ^ 0. 
The required matrix inversion, with a conjugate- gradient algorithm without any 
preconditioning, and the Fourier transform consume in average less than half the 
computing time of the Landau gauge fixing. 

4 Hypercubic artifacts 

The ghost propagator G(p) is a scalar invariant on the lattice which means that 
it is invariant along the orbit O (p) generated by the action of the isometry group 
if (4) of hypercubic lattices on the discrete momentum p = x (m, n 2 , n 3 , n 4 ) 
where the n^'s are integers, L is the lattice size and a the lattice spacing. The 
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general structure of polynomials invariant under a finite group is known from 
group-invariant theory. Indeed it can be shown that any polynomial function of 
p which is invariant under the action of H(4) is a polynomial function of the 4 
invariants p' n ' = a n ^uP^> n = 2,A,Q,8 which index the set of orbits. 

Our analysis program uses these 4 invariants to average the ghost propagator 
over the orbits of H{A) to increase the statistical accuracy: 

a*G L { p WJ±\jP\pW) = — d (P) (25) 

11 [P) " P eO(p) 

where ||0(p)|| is the cardinal number of the orbit 0(p). By the same token, one 
should always take the following real source 

S a p (x) =5 ab cos(p-x) (26) 
peo(p) 

rather than a single complex Fourier mode for studies of the ghost propagator 
in the infrared region. Indeed, after averaging over the gauge configurations and 
use of the translational invariance, one gets 

<Kim^))= £ j2((M- 1 )Z) e ~ tp, ' x+ip ' y 

p,p'eO(p) x,y 

= V\\0{p)\\a 2 G L {pWJ A \p^,pW) (27) 
By analogy with the free lattice propagator 

Go{p) ^ZM = ^( 1 + ^S + '")' where ft = r m (ir) (28 > 

it is natural to make the hypothesis that the lattice ghost propagator is a smooth 
function of the discrete invariants near the continuum limit, when ap^ 1 , V/i, 

G L (P 12] ,P 14 \P [6 \P 18] ) « G L ( P M, 0,0,0) +^^(^,0,0,0) + • • • (29) 

and Gl(p^ 2 \ 0, 0, 0) is nothing but the propagator of the continuuum in a finite 
volume, up to lattice artifacts which do not break 0(4) invariance. When several 
orbits exist with the same p 2 , the simplest method to reduce the hypercubic 
artifacts is to extrapolate the lattice data towards Gl(jP\ 0, 0, 0) by making a 
linear regression at fixed p 2 with respect to the invariant p' 4 ' since the other 
invariants are of higher order in the lattice spacing. The range of validity of 
this linear approximation can be checked a posteriori from the smoothness of the 
extrapolated data with respect to p 2 . 
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Choosing the variables p^ appropriate to the parametrization of a lattice prop- 
agator with periodic boundary conditions provides an independent check of the 
extrapolation. Indeed we can write as well 

G L (pM , , pM , P M ) = (pM , ^ (30) 

with the new invariants, again hierachically suppressed with respect to the lattice 
spacing, 



i>-' ■ Yl % ( 31 ) 



Gl and Gl are two different parametrizations of the same lattice data, but near 
the continuum limit one must also have 

G L (pM , pM , ^ , p^ ) « G L (pM , 0, 0, 0) + gf (^ , 0, 0, 0) + • • • (32) 

where G L (p [2 \ 0, 0, 0) and G L @ 2 \0, 0, 0) are the same function, the propagator 
of the continuum , of a different variable (again up to lattice artifacts which do 
not break 0(4) invariance). 

If one wants to include in the data analysis the points with a single orbit 
at fixed p 2 , one must interpolate the slopes extracted from Eqs (J29j) or (J32*Jl . 
This interpolation can be done either numerically or by assuming a functional 
dependence of the slope with respect to p 2 based on dimensional arguments. The 
simplest ansatz is to assume that the slope has the same leading behavior as for 
a free lattice propagator: 

p (p,2 ' -°' 0) = ^7 (ci + C2pI21) (33) 

The inclusion of 0(4)-invariant lattice spacing corrections is required to get fits 
with a reasonable x 2 . The quality of such two-parameter fits to the slopes, and 
the extension of the fitting window in p 2 , supplies still another independent check 
of the validity of the extrapolations. 

We have used Eqs. (J2HJ) and (j22J) to extrapolate our lattice data towards the 
continuum and determined the range of validity in p 2 of the extrapolations from 
the consistency of the different checks within our statistical errors. The errors on 
the extrapolated points have been computed with the jackknife method. TablesEl 
and E] summarize the cuts that have been applied to the data for the estimation 
of the systematic errors in the analysis of the next section. We have repeated the 
analysis of the gluon propagator j3] to study the sensitivity of the results with 
respect to the window in p 2 which has been enlarged considerably in our new 
data. The cuts for the lattice ghost propagator are stronger than for the gluon 
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p 


V 


Npoints 


" Pmin 


" Pmax 


x 2 


6.0 


16 4 


> 10 


< 1.30 


< 1.82 


< 1.4 


6.2 


24 4 


> 12 


< 1.30 


< 1.82 


< 1.1 


6.4 


32 4 


> 20 


< 1.40 


< 1.82 


< 1.3 



Table 3: Cuts on the lattice data for the gluon propagator. [ap m i n , ap max ] is 
the momentum window of a fit in lattice units and N po i nts is the number of data 
points in that window. 



p 


V 


Npoints 


& Pmin 


O Pmax 


x 2 


6.0 


16 4 


> 10 


< 1.30 


< 1.57 


< 1.0 


6.2 


24 4 


> 20 


< 1.30 


< 1.57 


< 1.0 


6.4 


32 4 


> 20 


< 1.00 


< 1.57 


< 1.0 



Table 4: Cuts on the lattice data for the ghost propagator. Columns have the 
same meaning as in TableEl 

lattice propagator because the statistical errors of the former are two to three 
times larger which make the continuum extrapolations less controllable. 

The number of distinct orbits at each p 2 increases with the lattice size and, 
eventually, a linear extrapolation limited to the single invariant breaks down. 
However there is a systematic way to include higher-order invariants and to extend 
the range of validity of the extrapolations. A much more detailed exposition of 
the controlling of systematic errors is in preparation, since our method has been 
largely ignored in the litterature where very empirical recipes are still in use. 



5 Data analysis 

The evolution equation of the renormalization constants of the gluon or ghost 
fields in a MOM scheme, with respect to the coupling constant h in an arbitrary 
scheme R (the index R is omitted but understood everywhere), can be written 
generically up to four-loop order: 

d\n{Z v ,MOM) _ _ , _ , 2 , — .3 , — 



cHn/i 2 



7^ + 71^ + 72^ + 73^ (34) 
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and the perturbative integration of Eq. (|34|) yields, to the same order, 



In — = log(/i) — + 



7o , , (A) 7i - Pi 7o) 



2 (/3q 2 7 2 -/3o/3i7i-(/3o/3 2 -/3i 2 )7o) 



2(31 (35) 
+ /> 3 (/3 3 7 3 - Pi Pi 7 2 + (Pa ft 2 - $ A) 7i 

+ (-/3 2 /5 3 + 2/5o/3 1 /3 2 -/3 1 3 )7o)-^ 
with the standard four-loop formula for the running coupling 



ft* ^ ft 2 i ft 4 * 2 27 /?? 4 

(36) 

, / M 2 
and t = log I — 

We now consider in turn the three renormalization schemes MS, MOMg and 
MOMc and fit the two parameters of Eqs. (|33|) and (f3Hj) to our extrapolated lattice 
data. FigureEl illustrates the typical quality of such fits. 



5.1 MS scheme 

The analysis in the MS scheme is summarized in Table El The statistical error is 
at the level of 1% for the gluon propagator and 2-3% for the ghost propagator, 
whereas the systematic error due to the extrapolations is around 3-5% and 5-10% 
respectively. The values of extracted from the gluon and the ghost propa- 
gators are consistent within these errors and within each order of perturbation 
theory. 

However, the three-loop and four-loop values, which are displayed in Table |H1 
with the physical units of Table El clearly confirm our previous result |3] that we 
are still far from asymptoticity in that scheme. 



5.2 MOMg scheme 

Table which summarizes the analysis in the MOMg scheme, shows that, at 
the lower /3's, we were not able to describe both lattice propagators at four-loop 
order with reasonable cuts and x 2 . This could be interpreted as an hint that 
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2 4 6 2 4 6 

p (GeV) p (GeV) 

Figure 2: Extrapolated lattice data at (3 — 6.4 for Z 3 (left) and Z 3 (right). The 
solid line is the fit at four-loop order in the MS scheme. The vertical dotted lines 
delimit the window of each fit. 



perturbation theory has some problems of convergence beyond three-loop order 
below 3-4 GeV. 

If we select the three-loop result as the best perturbative estimate of Aj^^j- 

and convert it to the MS scheme with the asymptotic one-loop formula, A^jg = 
0.346 Aj^jjj- , then we get the physical values quoted in Table |H1 which agree 
completely with previous values j3j. 

5.3 MOMc scheme 

The results of the analysis in the MOMc scheme are displayed in Table 01 We 

(3) 

still find that the three-loop and four-loop values of A J — are very much the 

L 1 MOMc J 

same both for the gluon propagator and for the ghost propagator. Thus the 
perturbative series seems again to become asymptotic at three-loop order in that 
scheme. 

Selecting the three-loop result as the best perturbative estimate of A^jg^ and 
converting it to the MS scheme with the asymptotic formula, A^g = 0.429 Aj^grj- , 
we get the physical values quoted in Table [TU1 
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p 


L 




A (2) 
MS, qluon 


X 


A (3) 

MS, qluon 


X 


A (4) 
MS,g«uon 


9 

X 


6.0 


16 


1.111 


0.336(3)1* 


1.3 


r 2 ^ — rr i 

0.265(3)1$ 


1.0 


0.225(2)1^ 


1.1 




24 


1.111 


0.332(3)1* 2 


0.6 


0.262(3)1^ 


0.5 


0.222(2)1* 


0.6 


6.2 


24 


0.907 


0.240(2)lg 


0.8 


0.185(2)1? 


0.8 


0.158(2)1| 


0.8 


6.4 


32 


0.760 


0.171(2)1^ 


1.4 


0.130(2)1}? 


1.4 


0.112(l)tjj 


1.4 


P 


L 




MS,qhost 


x 2 


MS,qhost 


x 2 


«A^ fc 

MS,qhost 


x 2 


6.0 


16 


1.039 


0.354(7)1^ 


0.5 


0.281(6)1^' 


0.5 


0.235(5)lf 


0.5 




24 


0.785 


0.325(6)1^ 


0.2 


0.259(5)1^ 


0.2 


0.217(4)1*3 


0.2 


6.2 


24 


0.693 


0.254(4)1$ 


0.4 


0.200(3)1^ 


0.4 


0.169(3)1^ 


0.4 


6.4 


32 


0.555 


0.193(2)1^ 


0.8 


0.150(2)1}* 


0.8 


0.128(2)1^ 


0.8 



Table 5: Fits of A^g from the gluon and ghost lattice propagators. The error 
in parenthesis is the statistical error corresponding to a window [ap m i n , ap max ] 
with the ap min quoted in the Table and the upper bound for ap max quoted in 
Tables El and E| respectively. 



(3 


A (3) 

MS. qluon 


MS. qhost 


Agi 

MS, qluon 


Ail 

MS,qhost 


6.0 
6.2 
6.4 


519(6)lf 
509(6)+^ 
476(7)1^ 


551(12)1^ 

550(8)11 

549(7)1^ 


441(4)1* 

435(6)1^ 

410(4)11 


461(10)1^ 

465(8)11 

468(7)12 



Table 6: Three-loop and four- loop physical values of A^jg in MeV extracted from 
Table El 



5.4 Scheme dependence 

The puzzling feature of Tables El El and [TQ] is the rather large dependence of the 
Aqcd scale upon the loop-order and the renormalisation scheme whereas, within 
any scheme, the values from the ghost and gluon propagators are rather consistent 
at each loop order and pretty independent of the lattice spacing. 

Let us consider again the evolution equation of the renormalisation constants 
of the gluon or ghost fields in a MOM scheme, with respect to the coupling 
in an arbitrary scheme R. We have 

d\n(Zr,MOM) _ v 1 d\n(Z rMO M) /0 ~n 
dln^ = = "2 d \nA R (3?) 

where Ar is the scale in scheme R. If we truncate the perturbative expansion at 
order n 



In 



(^f^)=c R ,oln(M + Ec«, fc 4 08) 
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R 
P 


T 




til V . — — 

MOMg,gluon 


V 2 
A 


aA (3) 

Cfci 1 . — — 

MOMg,.giuon 


V 2 
A 


til V — - — - 

MOMg,.^won 


V 2 
A 


6.0 


16 


1.039 


0.551(3)1| 


1.0 


0.477(3)lg 


1.2 








24 


1.014 


0.536(4)1- 


0.9 


0.464(3)11" 


0.9 






6.z 


O A 


0.693 


0.396(zJ_ 12 


1.0 


0.336(2)_ 15 


0.9 






6.4 


32 


0.555 


0.292(1)+^ 


1.3 


0.246(1)1^ 


1.4 


0.253(3)+^ 


1.6 




L 




MOMg,q/iosi 


x 2 


aA^ 

MOMf*„qhost 


x 2 


aA (4 ^ — 

MOMg,.ghost 


x 2 


6.0 


16 


1.039 


0.660(40)1^ 


0.4 


0.475(12)1^ 


0.5 








24 


1.014 


0.559(22)1^ 


0.2 


0.438(12)1- 


0.2 


0.408(17)1^ 


0.9 


6.2 


24 


0.693 


0.455(11)1^ 7 


0.3 


0.342(5)1^ 


0.6 


0.348(8)1$ 


1.0 


6.4 


32 


0.555 


0.333(4)1* 


1.2 


0.261(3)!^ 


0.9 


0.279(7)1- 


0.8 



Table 7: Fits of A^^^ from the gluon and ghost lattice propagators. The 
notations are the same as in Table 03 



(3 


WlS,qluon 


A®. 

MS,qhost 


6.0 
6.2 
6.4 


324(2)1^ 

320(2)1? 4 

312(1)155 


322(8)1^ 
326(5)+i 
331(4)1^ 



Table 8: Three-loop physical values of A^ in MeV extracted from Tabled 



the change in the effective scale A R \ or equivalently, the change in the coupling 
h R , induced by adding the contribution at order n + 1 is typically 

AA? ^ c R , n h R 



A W 2 lR {h 1 



(39) 



Now the dependence of the effective scale A R upon the coupling h R is given up 
to order 4 

2 In A< 4 > - In 1 ^InfBh ) ^ ~ ® h ^ " 2/Wa + ® h* 
2\nA R -\n» - — --] ll {p h R )-—^—h R h R 

(40) 

and, denoting the coefficient of order h R ~ 2 in that equation by —p Rn -i, the 
effective scales which describe a same coupling at order n and n + 1 are related 
by 

In-^-^n-i/^ 1 (41) 

A R 

Combining Eqs. (J39|) and (|41j) gives the relation between the effective scales which 
describe the renormalisation constants of the gluon or ghost fields in a MOM 
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R 
P 


j 




aA (2) 

Cii V - 

MOMc.giuon 


V 2 
A 


aA (3) 

(Xi V ■ 

MOMc, gluon 


V 2 
A 


aA (4) 

MOMc,gluon 


V 2 
A 


6.0 


16 


1.178 


0.482(6) 


1.3 


0.408(3)^4 


1.0 








24 


1.111 


0.468(5)^5 


0.5 


0.394(3)+6 i 


0.8 


0.411(7)12 


1.1 


o.z 


O A 
24 


o.yo< 


0.345(3)_ 10 


0.9 


U.2oo(Zj_ 6 


0.9 


0.2y2(oJ_ 5 


O.O 


6.4 


32 


0.589 


0.255(1)1^ 


1.5 


0.205(1)+^ 


1.6 


0.212(1)+^ 


1.6 




L 




aA^ 

MOMcqhost 


x 2 


aA^ 

MOMc, qhost 


x 2 


MOMc,ghost 


x 2 


6.0 


16 


0.962 


0.489(6)1^ 


0.8 


0.437(11)1^ 


0.4 








24 


1.047 


0.459(6)1^ 


0.5 


0.408(8)1!? 


0.5 


0.398(14)1^ 


0.3 


6.2 


24 


0.740 


0.367(7)+^ 


0.4 


0.308(7) + ^ 6 


0.2 


0.303(9)1( 4 


0.2 


6.4 


32 


0.589 


0.280(5)+i 


0.6 


0.225(5)1- 


0.6 


0.224(5)1^ 


0.6 



Table 9: Fits of Ajjj^r from the gluon and ghost lattice propagators. The nota- 
tions are the same as in Table El 



p 


A^ 

WS,qluon 


A (3) 

MS, ghost 


6.0 
6.2 
6.4 


345(3)11 
341(2)1? 
323(2)1^ 


369(9)1°; 

364(8)1^ 

354(8)1^ 



Table 10: Three-loop physical values of Aj^g in MeV extracted from Table H3 



scheme at order n and n + 1 

^ r = exp- 5 ^,„_ 1 + — j fts (42) 

Figure El displays the behavior of this ratio for the gluon and ghost propaga- 
tors in the three schemes as a function of the momentum p for n = 2 and n = 3. 
The couplings are taken from the fits at /? — 6.4. There is a pretty good quali- 
tative agreement with Tables IBl IH1 and ITU1 which confirms the overall consistency 
with perturbation theory of the lattice data for the gluon and ghost propagators 
within any renormalization scheme. 

The scheme dependence of the Aqcd scale can also be analyzed with Eq. (|4TIjl : 

< ) - exp (,L^_ I LV4^ + '''} (43) 



y \2Po\h Rl h R J 2(31 h 



Figure |U shows the behavior of the ratios ( „y s and ^° Mc , as a function of the 

MOMg MOMg 

momentum p at each order of perturbation theory. The couplings are taken from 
the fits of the gluon propagator at (3 — 6.4. 

Clearly, the limiting values of these ratios are not the asymptotic values. If we 
replace in Eq. (|4"3|) the coupling Jir 2 by its perturbative expansion with respect 
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to h Rl 



h R2 =h Rl +Y,r k h% 1 (44) 



k=l 



then the ratios do of course tend towards the asymptotic values exp | ^-J- I . The 

disagreement with respect to the perturbative expansion is not a problem with 
the lattice data or with the numerical analysis. Indeed the fits do a very good job 
at extracting a well-behave coupling as illustrated in Fig. |5] which displays the 
dimensionless scales aA^L aA^ — - and aA^ — as a function of the momentum 

MS' MOMg MOMc 

p, using Eq. ()4U|1 with the fitted couplings at (3 = 6.4 from the ghost and gluon 
propagators. Z , the other fitted parameter of Eq. is nearly independent, 
within a few percent, of the renormalisation scheme as it should in the absence 
of truncations. It follows that the difficulty to reproduce the asymptotic ratios 
between the scales of different renormalization schemes, is mainly a consequence 
of the truncation of the perturbative series of the renormalization constants of 
the gluon and ghost propagators. 

We can substantiate this claim, and estimate the rate of convergence, by the 
following exercise. We solve h R2 in terms of h Rl using Eq. (|38j) at four-loop order 

In ( T '™ OM j = c R2i0 ln(/iK 2 ) + y2c R2tk h k R2 = c Ru0 \n(h Rl ) + y^c Rl , k h k Ri (45) 
\ Z ° J k=i ' k=i 

Then we plug the solution into Eq. (|43jl. Figure |H1 shows the behavior of the 

aIIL 

corresponding ratios, { ™ s and ^° Mc , as a function of the coupling h^g and 

MOMg MOMg 

^momc res P ec tively. The effect of the truncation of the perturbative series is 
manifest for the MS scheme and gives the right order of magnitude of what is 
actually measured in Tables El and 



6 Conclusion 

We have shown that the lattice formulation of the ghost propagator has the ex- 
pected perturbative behavior up to four-loop order from 2 GeV to 6 GeV. We 
have been able to go beyond the qualitative level and to produce quantitative 
results for the scale A^ which are pretty consistent with the values extracted 
from the lattice gluon propagator. We have understood the strong dependence 
of the effective Aj^g scale upon the order of perturbation theory and upon the 
renormalisation scheme used for the parametrisation of the data. The perturba- 
tive series of the MOM schemes seem to be asymptotic at three-loop order in the 
energy range we have probed whereas the MS scheme converges very slowly. If we 
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assume that all perturbative series remain well behaved beyond four-loop above 
4 GeV, then we get A^jg & 320 MeV with a 10% systematic uncertainty. The 
statistical errors are at the 1% level. This value is also in pretty good agreement 
with the values of A^jg extracted from the three-gluon vertex in a MOM scheme 
at three- loop order [2T], at the same /3's and with the same lattice sizes. On the 
other hand it exceeds by 20% the value obtained from the same vertex at /3 — 6.8 
on a 24 4 lattice. This discrepancy motivated the introduction of power correc- 
tions which are successful in describing the combined data of the three-gluon 
vertex j!4j . We will show in a forthcoming paper how the power corrections can 
be unraveled from the lattice propagators alone. 

The value quoted above exceeds also by about 30% the previous determina- 
tions of the QCD scale in the quenched approximation based on gauge-invariant 
definitions of the strong coupling constant [221 (take note, for comparison 
purposes, that our physical unit corresponds to the force parameter r |23j set 
approximately to 0.53 fm). However there is also an uncertainty due to the use 
of the asymptotic one-loop relation between A^jg and the A L 's. For illustration, 
let us consider the determination of Aj^g using lattice perturbation theory up to 
three-loop order with the Wilson action ■ It is possible to estimate the rate of 

A (3) 

convergence of the ratio — ^j- as a function of the bare lattice coupling h L = 

by inserting the perturbative expansion of %jg into Eq. (|33J) . Figure U\ displays 

A (3) 

the evolution of this ratio and also of the ratio — Rv- for the so-called "boosted" 

MS 

lattice scheme which re-express the lattice perturbative series as a function of the 
coupling hu — hij/ (plaq). The mere truncation of the perturbative series intro- 
duces an uncertainty on the absolute scale of the lattice schemes which could be 
as large as 30% in the range of (3 studied in these simulations. 

No strategy can fix the scale A QCD to an accuracy better than the uncertainty 
entailed by the truncation of the perturbative series in the conversion to the 
MS scheme. We have shown that this error can be larger than the main well- 
known sources of systematic errors which come from setting the scale a -1 and 
from the continuum extrapolation. If we aim at reducing below 10% the error in 

the conversion of the MOM schemes to the MS scheme, then a look at Figure |H1 
shows that we need to apply a cut at 6 GeV. Such an analysis would require 
simulations at (3 = 6.6 and (3 = 6.8 on 48 4 and 64 4 lattices respectively, to work 
at fixed volume and minimize finite-size effects. The existence of several lattice 
observables, gluon propagator, ghost propagator, three-gluon vertex, from which 
one can extract independent values of the scale Aqcd, an advantage of the Green 
function approach, should then allow to disentangle unambiguously the effects of 
the truncation of the perturbative series from the non-perturbative corrections, 
and to get a value of Aj^g at a true 10% accuracy. 
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2 4 6 23456 

p (GeV) p (GeV) 



Figure 3: fl (n) for n = 2 (dashed lines) and n = 3 (solid lines), for the gluon 

propagator in the MS scheme (a), MOMg scheme (c) and MOMc scheme (e), and 

for the ghost propagator in the MS scheme (b), MOMg scheme (d) and MOMc 
scheme (f). 
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Figure 4: A( „y s (left) and A ^ff Mc (right) for n = 2, n = 3 and n = 4. The solid 

MOMg MOMg 

lines are the plots of Eq. ()43j) with the fitted couplings whereas the dashed lines 
are the plots with Eq. (|44|) . Horizontal lines are the asymptotic values. 
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Figure 5: aA^^L aA^^— and aA^ — from the gluon propagator (solid lines) 

& MS' MOMg MOMc & f f o \ ) 

and from the ghost propagator (dashed lines) at (5 — 6.4, as a function of the 
momentum through Eq. ([40)) . 



24 




Figure 6: { ™ s (solid line) as a function of hy^ and ^ ) OMc (dashed line) as a 

MOMg MOMg 

function of h^~r^ c . The vertical lines delimit the values spanned by (dashed) 
and h^z-r (dotted) in the fits of the gluon propagator at (3 — 6.4. 
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Figure 7: -^y (lower solid line) as a function of h L and -By- (upper solid line) as 

MS" MS 

a function of ha- The vertical lines (dotted) delimit the values spanned by 
and hu in the simulations of |2S] (5.7 < (3 < 6.9). The dashed horizontal lines 
are the asymptotic values. 
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